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Abstract: 

The spontaneous symmetry breaking in the Ising model can be equivalently described 
in terms of percolation of Wolff clusters. In 0(n) spin models similar clusters can be 
built in a general way, and they are currently used to update these systems in Monte 
Carlo simulations. We show that for 3-dimensional 0(2), 0(3) and 0(4) such clusters are 
indeed the physical 'islands' of the systems, i.e., they percolate at the physical threshold 
and the percolation exponents are in the universality class of the corresponding model. 
For 0(2) and 0(3) the result is proven analytically, for 0(4) we derived it by numerical 
simulations. 



1 Introduction 

The possibility to interpret the critical behaviour of dynamical systems in terms of perco- 
lation of geometrical structures of the system has always had a great appeal in the study 
of critical phenomena JH - Attempts in this direction were already done in the 70 's, 
when one began to study the behaviour of clusters of nearest-neighbour like-signed spins 
in the Ising model. It turned out that in two dimensions these elementary site percolation 
clusters indeed undergo a geometrical transition exactly at the critical threshold of the 
Ising model 0. This result, which is not valid in three dimensions [Q, is anyway not 
so appealing, because the critical exponents derived by the percolation variables do not 
coincide with the Ising ones 0. The correspondence between the geometrical and the 
thermal phenomenon is therefore only partial. 

The problem was solved by A. Coniglio and W. Klein , making use of a different def- 
inition for the clusters. Such definition had already been used by Fortuin and Kasteleyn 
to show that the partition function of the Ising model can be rewritten in purely geomet- 
rical terms as a sum over clusters configurations 0. According to the Fortuin-Kasteleyn 
prescription, two nearest-neighbouring spins of the same sign belong to the same cluster 
with a probability p = 1 — exp{—2(3) {j3 = J/kT, J is the Ising coupling). Coniglio and 



Klein showed that the geometrical transition of these clusters leads to the required 
critical indices of the Ising model, threshold and exponents. 

The clusters of the Monte Carlo cluster update introduced by U. Wolff [0 for 0(n) spin 
models coincide with the Fortuin-Kasteleyn ones when n=l (which is just the Ising model). 
The 0(n) models without external field in three space dimensions {n>2) undergo a phase 
transition due to the spontaneous breaking of the continuous rotational symmetry of 
their Hamiltonian. Such models are very interesting: some physical systems in condensed 
matter physics are directly associated to them. The three-dimensional 0(3) model is the 
low-temperature effective model for a bidimensional quantum antiferromagnet The 
0(2) model in three dimensions is known to be in the same universality class as superfluid 
^ife. 0(n) models are also very useful to study relativistic field theories. The 0(4) model 
in three dimensions has been conjectured to be in the same universality class as the 
finite-temperature chiral phase transition of QCD with two flavours massless quarks 0. 

The general definition of Wolff clusters for 0(n) spin models inspired this work. Can 
one describe the critical behaviour of 0(n) without field in three dimensions in terms of 
the percolation of these clusters, like in the case n=l? We will show that this is indeed 
true at least for 0(2), 0(3) and 0(4). The fact that the Wolff clusters percolate at 
the physical critical point was recently proven analytically for 0(2) and 0(3) |T0|, |lT|. 
Although nothing about the exponents was mentioned, we will show that starting from 
some relations established in 0, [ll[] it is also possible to deduce the equality of the critical 
exponents for 0(2) and 0(3). We have also performed computer simulations on 0(2) and 
0(4) in order to illustrate this result for 0(2), and to prove it numerically for 0(4). 

2 0(n) models and Wolff clusters 

The 0(n) spin models with no external magnetic field have the following Hamiltonian: 



H = -jY,s,Sy (2.1) 

(iJ> 

where i and j are nearest-neighbour sites on a d-dimensional hypercubic lattice, and Si is 
an n-component unit vector at site i (J is the coupling). The partition function of these 
models at the temperature T is 



Z(T) = j I?[s]exp{/3 5^SiSj} (2.2) 



(i,j> 

where (3 = J/kT and the integral is extended over all spin configurations {s} of the 
system. In three dimensions the 0(n) models undergo a second-order phase transition. 
The order parameter of this transition is the normalized magnetization M = ^ Yli (V 
is the lattice volume). 

Numerical simulations of 0(n) models became much quicker and more effective after 
U. Wolff introduced a Monte Carlo algorithm based on simultaneous updates of large 



clusters of spins, generalizing the Swendsen Wang algorithm |12] to the continuous-spin 



case. This algorithm has the remarkable advantage that it eliminates the problem of 
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critical slowing down, an effect that makes simulations around criticality very lengthy 
with traditional local methods (Metropolis, heat bath). The Wolff algorithm can be 
basically divided in two phases: 

1) a cluster of spins is selected; 

2) the spins of this cluster are "flipped", i.e. they are reflected with respect to some 
defined hyperplane. 

For details of the flipping procedure, see 0. Here we are interested in the way to 
build up the clusters. We can split this procedure in two steps: 

a) choose a random n-component unit vector r; 

b) bind together pairs of nearest-neighbouring sites i,j with the probability 

p{i,j) = 1 — exp{min[0, — 2/3(si-r)(sj-r)]}. (2.3) 

From this prescription it follows that if the two spins at two nearest-neighbouring sites i 
and j are such that their projections onto the random vector r are of opposite signs, they 
will never belong to the same cluster {p{i,j) = 0). The random vector r, therefore, divides 
the spin space in two hemispheres, separating the spins which have a positive projection 
onto it from the ones which have a negative projection. The Wolff clusters are made out 
of spins which all lie either in the one or in the other hemisphere. In this respect, we can 
again speak of 'up' and 'down' spins, like for the Ising model. In addition to that, the 
bond probability is local, since it depends explicitly on the spin vectors Si and sj, and not 
only on the temperature like the Fortuin-Kasteleyn factor. 

The analogies with the Ising model are however clear, motivating the attempt to study 
the percolation properties of these clusters. 



3 Percolation exponents for 0(2) and 0(3) 

In |Tl|, the random cluster representations of 0{n) models, n = 2,3 have been 
derived (and exploited) through the Fortuin-Kasteleyn transformation of the Hamilto- 
nians and similar results were obtained in |T^ for the continuous (or classical) spin model 



[p!5| . Wolff random cluster probability distributions for these models have been stud- 
ied and several monotonicity properties of these distributions (FKG properties) [|T^ have 
been established leading to the proof of the equivalence between the onset of magnetic 
ordering in the 0{n) models, n = 2, 3 and percolation in the corresponding random Wolff 
cluster models. From the results stated above follows in a natural way the equality of the 
critical thermal and geometrical exponents. 

In what follows we focus our attention on two variables: 

• The percolation strength P, defined as the probability that a lattice site picked 
up at random belongs to the percolating cluster. P is the order parameter of the 
percolation transition. 

• The average cluster size S, defined as 

5 = 5f^- (3.1) 

Here, rig is the number of clusters of size s per lattice site (i.e. divided by the lattice 
volume) and the sums exclude the percolating cluster (see 
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We stress that in some cases one speaks of average cluster size referring to the average 
size of all clusters, including the percolating one. This definition makes the variable 
infinite above the percolation threshold. Below the threshold the two definitions obviously 
coincide and therefore they share the same critical exponent. 

In the references |TT], |T3|, |T^, it was proved that, if P{T) is the percolation strength 
in the random Wolff cluster representation and m{T) is the magnetization in the 0{n) 
spin system, then there exists a function c(T) G C°°(JR'^) such that 

P(T) > m(T) > c(T) P(T) (3.1) 

Similarly, denoting S{T) the average size of the Wolff clusters and x{T) the (linear re- 
sponse) susceptibility, it was also proved that 

SiT) > xiT) > c(T)2S(T) (3.2) 

where c(T) is the same C°°(1R) function as in (p.l|) . 

From scaling theory (see ll^), near criticality, the susceptibility is believed to behave 
according to the following law 

X(T) ^ (T-T,)-^ (3.3) 

and the average mean cluster size of Wolff clusters should follow the law 

S{T) ^ \p{T)-p{n)r' (3.4) 

where p{T) is the bond occupation probability in the Wolff random cluster model (i.e. 
the Coniglio-Klein bond probability) given by p{T) ~ 1 — exp{—a/T), where a does 
not depend on T. 

From [|T0], |Tl], |T^, |1^] the equality of the critical exponents 7 and 7' follows readily . 



Indeed, because of monotonicity, taking the logarithm in ( p.3[ ) and (|3.4| ), and using ( p.2|) 
one gets 

-7' log I p(T) - p(T,) I > -7 log(T - T,) > 2 log c(T) - 7' log | p(T) - p{T,) \ (3.5) 
which reduces to 

y>-, >y !i2l£(I)_ (3.6) 

log I p{T) - p{Tc) I log I p{T) - p{Tc) I 

Now when T and since c(T) G C°°(1R), the last term vanishes and it is easy to see 

that log(T — Tc)/log I p{T) — p{Tc) \ 1. When T T+ we get the result 7 = 7'. 

Using (|3.1|) and the scaling behaviours of the magnetization and the percolation 
strength in terms of their critical exponents (3 and (]' respectively, one can show (fol- 
lowing the same lines as before) that 13 = 13' when T = p{T) p{Tc)^, where p{T) 
is again the Coniglio-Klein bond probability already defined above. 

Namely, the percolation strength is believed to behave |18[ as a function of the ele- 
mentary bond occupation probability p according to the following law 



P{p) ^{p-PcY P^P{T) (3.7) 



4 



whereas, the magnetization should behave as 



m(T) ~ (T,-T)^ (3.8) 

then using the same procedure as below, we are led to the following expression 
fl'<fl log(^c - T) \ogc{T) 

^ - ^ iog(p(r) - p(T,)) - ^ iog(p(T) - Pin)) ^ • ^ 

which, using logiT^-T)/ log(p(T) -piT^)) 1 and log c{T)/ log {p{T) -piT^)) 0, 
gives 

13 = 13'. 

as claimed. 

About eventual extensions of analytical results to other spin models, we mention that 
the equivalence of the thermal and the geometrical phase transitions is true for Z{3) and 
Z(4) models as a consequence of the equivalences already established for the continuous 
spin Ising model and 0(2) |TU|. As for 0{n) models with n > 3, we remark that the 



proof for 0(3) was derived starting from the result for 0(2). However, some additional 
conditions remain to be proved in order to be able to extend the result from 0(?7. — 1) to 
0{n) if n > 3. 



4 Numerical Results 

We have investigated numerically the 3-dimensional 0(2) and 0(4) models performing 
computer simulations for several lattice sizes. The Monte Carlo update was performed by 
the Wolff algorithm, described in Section 2. At the end of an iteration, the percolation 
strength P and the average cluster size 5* were measured. This has been done for each of 
the models using two different approaches. 

The first approach is the traditional one, based on a complete analysis of the lattice 
configuration. Once we have the configuration we want to analyze, we build Wolff clusters 
until all spins are set into clusters. We assign to P the value zero if there is no percolat- 
ing cluster, the ratio between the size of the percolating cluster and the lattice volume 
otherwise. We calculate S using the standard formula ( |3.1| ). The operative definition of 
percolating cluster was taken as follows. We say that a cluster percolates if it spans the 
lattice from a face to the opposite one in each of the three directions x, y, z. We made 
this choice to reduce the possibility that, due to the finite size of the lattices, one could 
find more than a spanning cluster making ambiguous the definition of our variables^. 

In this approach we have used free boundary conditions. 

The second approach is based on a single-cluster analysis. Basically one studies the 
percolation properties of the cluster built during the update procedure. For the cluster 
building we have considered periodic boundary conditions. Suppose that Sc is the size of 

^In three dimensions even this definition of spanning cluster does not exclude the possibility of having 
more than one of such clusters for the same configuration. Nevertheless the occurrence of such cases is 
so small that we can safely ignore them. 
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the cluster we built. If it percolates, we assign value one to the strength P and zero to 
the size S\ otherwise, we write zero for P and s^. for S. These definitions of P and S 
look different from the standard definitions we have introduced above, but it is easy to 
see that they are instead equivalent to them. 

In fact, we build the cluster starting from a lattice site taken at random. In this 
way, the probability that the cluster percolates (expressed by the new P) coincides with 
the probability that a site taken at random belongs to the percolating cluster (standard 
definition of P). As far as the average cluster size is concerned, we can repeat the same 
reasoning: the probability that the cluster we built is a non-percolating cluster of size Sc 
is just the probability that a randomly taken lattice site belongs to a non-percolating 
cluster of size Sc] Wg^ is given by 

Ws, = Sc . (4.1) 

Because of that, whenever we get a non-vanishing size Sc, such value will be weighted by 
the probability Wg^ in the final average S, which is then given by the following formula: 

S = ^Ws.Sc = ^ris^Sc^, (4.2) 

where the sum runs over the non-percolating clusters. We notice that Eq. ( [4.2|) coincides 
with Eq. (|3.1|) , apart from the denominator ^^^^s s, which is just the density of the sites 
belonging to finite clusters. Since this term does not contribute to the divergence of the 
average cluster size, the power law behaviour of the two 5"s at criticality is identical, so 
that the critical exponent 7 is the same in both cases. 

As we have said, in the second approach we select a single cluster at a time from the 
whole configuration. Because of that we have now some freedom of choosing the definition 
of percolating cluster, as we do not risk, like in the first approach, to find more spanning 
structures. We say that the cluster percolates if it connects at least one face with the 
opposite one. 

In this way, also the definitions of percolating clusters are different in the two ap- 
proaches. This certainly influences the results on finite lattices, but has no effects on the 
infinite- volume properties we are interested in. In fact, it is known that one can have 
at most a unique spanning cluster above the critical density Pc (in our case below the 
critical temperature Tc). Exactly at Pc (T^) there is a finite probability to have more than 
a spanning cluster [|T^. So, the two different definitions of percolating cluster we have 
adopted can lead to differences between the infinite-volume values only at the critical 
point Pc (Tc). But the critical exponents are, of course, not influenced by that, as they 
are determined by the behaviour of the percolation variables near the critical point, not 
exactly at Pc (Tc). 

The second approach has the advantage that it does not require a procedure to reduce 
the configuration of the system to a set of clusters; on the other hand, since it gets the 
information out of a single cluster, it requires a higher number of samples in order to 
measure the percolation variables with the same accuracy of the first method. Neverthe- 
less, the iterations are faster due to the simpler measurement of observables, and are less 
correlated than in the first approach, since only a (random) limited region of the lattice 
is considered in each sample. We find that both methods are efficient, and that it is 
important to be able to compare results obtained in two such different ways. 
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For our numerical investigations we have also made use of another variable which can 
be extracted from the percolation strength P. On a finite lattice there is at any temper- 
ature (3 a well defined probability of having a spanning cluster. We call it percolation 
cumulant and indicate it 7,.. When the size of the lattice goes to infinity, 7^ as a function 
of P approaches a step function: it is always zero below (3^ and always one above it. To 
get the finite-size curves out of our measurements we must basically see how often we 
found a percolating cluster (Pj^O) for a definite lattice size and a temperature p. 

To evaluate the thresholds and the exponents we have adopted finite-size-scaling tech- 
niques. We consider the general finite-size-scaling prediction for an observable O 

Oit,L) = L^/^go(Li/^t) , (4.3) 

where t = T — Tc , L is the linear dimension of the lattice, Qo is a universal function 
and the exponent p is related to the critical behavior of O at infinite volume. Following 
the definitions given in Section 3, we have p = 7' for the observable S and p = —(5' 



for P. For the percolation cumulant 7^ we have p = |18|, which means that 7,. curves 
corresponding to different lattice sizes cross at the critical point: for t = and p = 0, in 
fact, the observable of Eq. ( [4.3|) is not L-dependent. 

Figs. 1 and 2 show 7^, curves for 0(2) and 0(4), respectively. The agreement with 
the physical thresholds (dashed lines) is clear. Moreover, we could already get indications 
about the class of critical exponents of our clusters. In fact, if one knows the critical point 
and the exponent u, a rescaling of 'jr as a function of {T — Tc)L^^'^ should give us the same 
function for each lattice size (see Eq. ^^) . Figs. 3 and 4 show the rescaled percolation 
cumulant curves for 0(2), using (3c = 0.45416 and two different values of the exponent z/, 
respectively the 0(2) value and the random percolation one. 
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Figure 1. Percolation cumulant as function of /? 
for 0(2) and five lattice sizes. The dashed line 
indicates the position of the thermal threshold 
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Figure 2. Percolation cumulant as function of P 
for 0(4) and six lattice sizes. The dashed line 
indicates the position of the thermal threshold 




The scaling we get in correspondence of the 0(2) value is remarkable. In Figs. 5 
and 6 we repeat the same analysis for 0(4) {Pc = 0.9359); also here it is clear that the 
percolation exponent u is in agreement with the 0(4) value. (We have considered values 
for the 0(2) and 0(4) models from Refs. 



2g and [|T], m respectively.) 
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To determine more precisely the critical point we have used the scaling relation ( [4 .31 ) 
for the variables 5* and P. By plotting O as a. function of L at the critical temperature, 
we can obtain the exponents' ratio p/u directly from the slope of the data points in a 
log- log plot. 

We concentrated ourselves on the critical regions that we localized through the perco- 
lation cumulant and performed more simulations for several f3 values looking for the /3's 
for which we get the best for the linear fit of the data points in a log- log plot. The 
results for 0(2) and 0(4) are in Tables I and II, respectively. 
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Figure 3. Rescaled percolation cumulant for 0(2) 
using Pc = 0.45416 and the 0(2) exponent i' — 
0.672. 
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using (3c = 0.9359 and the 0(4) exponent — 
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f3/u 


7/z/ 


Percolation results 


0.45418(2) 


0.516(5) 


1.971(15) 


Thermal results 20| 


0.454165(4) 


0.5189(3) 


1.9619(5) 



Table I. Comparison of the thermal and percolation thresholds and exponents for 0(2) . 







f3/u 


7/z/ 


Percolation results 


0.93595(3) 


0.515(5) 


1.961(15) 


Thermal results 


0.93590(5) |21[ 


0.5129(11) g 


1.9746(38)0 



Table II. Comparison of the thermal and percolation thresholds and exponents for 0(4). 



The agreement with the physical values in Refs. |2^, 0, ^ is good. So far we have 
presented the results obtained using the first approach. The results derived using the 
second approach are essentially the same; besides, we observe an improved quality of 
the scaling, mainly because of the use of periodic boundary conditions, which reduce 
considerably the finite-size effects. In particular we show in Figs. 7, 8 the scaling of S 
and P at the thermal thresholds reported in Refs. pO|, We observe very small finite- 
size effects (lattices of L > 20 are used in the fits), especially for the 0(2) case, which is 
in contrast to what is observed for thermal observables ||23|. The slopes of the straight 
lines are in agreement with the values of the thermal exponents' ratios /5/z/, 7/z/. 

We thus confirm numerically the equivalence found in Section 3 for the 0(2) case, and 
verify that it holds also in the 0(4) case. 
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Figure 7: Finite-size-scaling plot at for the percolation observable P as a function 
of the lattice size L. The slopes in the plots correspond to (3'/v' = 0.521(3), 0.513(6) 
respectively for 0(2) and 0(4). Error bars are one standard deviation. 
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S for the 0(2) Model 
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Figure 8: Finite-size-scaling plot at Tc for the percolation observable 5" as a function of the 

lattice size L. The slopes in the plots correspond to 7'///' = 1.97(1), 1.99(1) respectively 
for 0(2) and 0(4). The two curves look surprisingly similar to each other. Error bars are 
one standard deviation. 

5 Conclusions 

In this work we have shown that the spontaneous breaking of the continuous rotational 
symmetry for the 3-dimensional 0(2), 0(3) and 0(4) spin models can be described as 
percolation of Wolff clusters. For 0(2) and 0(3) the result was proven analytically, for 
0(4) it was derived by means of lattice Monte Carlo simulations. In all cases, the number 
n of components of the spin vectors s does not seem to play a role; the result is thus likely 
to be valid for any 0(n) model. 
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